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PART  I.     ANALYSIS  OF  A  MULTI-STAGE  AXIAL  COMPRESSOR 
WITH  HIGH  REACTION  BLADING 


ABSTPvACT  FOR  PART  I 

An  analysis  of  a  multi-sta^e  axial  compressor  with  high  reac- 
tion blading  is  carried  out.     The  usual  methods  of  linearizing  the  equa- 
tions for  flow  through  the  compressor  fail  for  this  type  of  blading.     A 
numerical     solution  of  the  non-linear  equations  is  worked  through  for 
off-design  operation.     Flow  velocities  through  the  first  six  stages  are 
calculated  for  the  off-design  flow  rates  in  order  to  insure  that  the  flow 
does  indeed  come  to  a  steady- state  repeating  pattern.     Entering  flow 
angles  for  both  rotor  and  stator  are  calculated  for  this  repeating  flow 
condition.     The  off-design  incidence  angles  for  the  stators  at  low  flow 
rates  indicate  a  possible  deterioration  of  efficiency. 


SUMMARY  FOR  PART  I 

The  pressure  rise  per  stage  in  an  axial  flow  compressor  is  pro- 
portional to  the  square  of  the  incidence  velocity  relative  to  the  blades. 
Almost  all  axial  compressors  for  gas  turbines  were  designed  with  the 
whirl  in  the  same  direction  as  the  rotor  blade  speed.     This  choice  of 
design  flow  gave  high  efficiencies  when  the  blades  were  designed  prop- 
erly.    Until  the  introduction  of  transonic  blading,   a  high  efficiency  was 
obtained  only  when  the  Mach  number  relative  to  the  blade  was  less  than 
unity  --  usually  in  the  neighborhood  of  0.  75.     Relative  velocities  of  this 
magnitude  did  not  require  whirl  opposite  to  the  direction  of  rotor  blade 
speed  and,   in  fact,   the  matching  with  a  turbine  was  easier  with  whirl  in 
the  same  direction  as  the  rotor  speed. 

An  axial  compressor  for  a  gas  of  low  molecular  weight,  for  in- 
stance helium,   can  be  designed  with  blading  similar  to  the  axial  compres- 
sor for  air.     However,   the  pressure   rise   per  stage  is  proportional  to 
the  square  of  the  Mach  number  of  the  relative  velocity.     The  rotor  blade 
rotational  speed  cannot  be  increased  because  of  strength  limitations,   and 
hence  the  square  of  the  relative  Mach  number  and  therefore  the  pressure 
ratio  are  decreased  by  a  factor  approximately  proportional  to  the  ratio 
of  molecular  weights.     Thus,   the  number  of  stages  required  for  a  com- 
pressor using  helium  would  be  approximately  seven  times  that  for  air 
for  the  same  overall  pressure.     Clearly,   higher  relative  velocities  are 
required  to  decrease  the  number  of  stages  for  a  gas  of  lower  molecular 
weight.     An  obvious  method  of  doing  this  is  to  introduce  the  whirl  in  the 
opposite  direction  to  rotor  blade  speed. 

A  vortex  flow  pattern  was  chosen  as  the  design  condition:    then 
the  flow  is  truly  two-dimensional;  that  is,   there  is  no  radial  shift  of  the 
streamlines  as  the  fluid  proceeds  downstream  through  the  compressor. 


This  considerably  simplifies  the  analysis  of  the  flow  at  the  design  flow 
rate.     Off-design  calculations  must  still  be  carried  out  in  order  to  in- 
vestigate the  compressor  performance  and  to  indicate  how  rapidly  losses 
might  increase  as  the  compressor  flow  rate  increases  or  decreases  from 
the  design  value.     For  conventional  multi-stage  compressors  with  rela- 
tively low  ideal  pressure  coefficients,   this  analysis  is  usually  carried 
out    by  linearizing  the  governing  equations  for  three-dimensional  flow 
through  the  compressor,   and  is  quite  satisfactory.     However,   with  whirl 
against  the  rotor  blades  and  its  correspondingly  higher  pressure  coef- 
ficient,  this  type  of  linearization  fails  completely.     It  is  necessary  to  al- 
low for  the  streamline  shift  to  avoid  singularities  in  the  solution.     Thus, 
the  full  set  of  non-linear  equations,   taking  into  account  the  radial  shift 
of  the  streamlines,   must  be  used.     In  order  to  calculate  the  repeating 
flow  pattern  characteristic  of  multi-stage  compressors,   it  was  necessary 
to  numerically  determine  the  flow  pattern  through  the  first  four  or  six 
stages. 

Although  detailed  calculations  of  efficiency  are  scarcely  possible, 
the  results  of  the  analysis  can  be  used  to  give  indications  of  probable 
losses.     The  degree  of  turning  required  by  the  blades  appears  to  be  with- 
in permissible  limits.     The  off-design  incidence  angles  for  the  rotors 
deviate  very  little  from  the  design  values.     The  stator  incidence  angles 
are  more  severe,  and  might  produce  a  loss  in  efficiency  over  part  of  the 
stator  blades. 


I.     INTRODUCTION 

Conventional  subsonic  axial-flow  compressors  designed  to  use 
air  as  a  working  fluid  are  not  necessarily  aiitable  for  gases  of  lower 
molecular  weight.     For  a  particular  closed  cycle  system,  where  a 
working  fluid  other  than  air  might  be  desirable,   the  number  of  stages 
required  will  vary  greatly  with  the  choice  of  gas.     The  number  of  stages 
is  directly  proportional  to  the  specific  heat  at  constant  pressure  and  is 
inversely  proportional  to  the  blade  speed  and  the  change  in  tangential 
velocity  of  the  gas  through  a  stage.     The  blade  speed  is  fixed  by  the 
strength  of  the  blade  material  and  is  generally  held  as  high  as  possible. 
With  present-day  air  compressors,  there  is  a  maximum  amount  of 
tangential  velocity  change  that  can  be  accomplished  efficiently.     Thus, 
for  a  typical  compressor  designed  to  use  air,  the  number  of  stages  re- 
quired would  be  increased  approximately  five  times  if  a  gas  such  as 
helium  were  used  and  the  same  temperature  ratio  across  the  compres- 
sor were  to  be  maintained.     This  report  describes  a  blade  design  that 
will  hopefully  reduce  this  high  number  of  stages  required  for  helium  to 
approximately  one-half  the  number  presently  used. 
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II.     THEORY 

Unlike  conventional  multi-stage  compressors,  the  design  of  the 
stator  blade  is  such  as  to  give  an  accelerating  flow  through  the  stator 
rows.     The  rotor  blades  then  turn  the  flow  back  to  the  axial  direction. 
This  type  of  design  has  been  used  in  some  single-stage  fans.     The  nota- 
tion to  be  used  in  this  report  is  illustrated  in  figure  1. 
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Tigure  1. 


The  following  list  of  symbols  will  be  used  throughout  this  report: 
c      =    axial  component  of  velocity 
c      =    tangential  component  of  velocity 
c      =    absolute  velocity 
w     =    velocity  relative  to  the  rotor 
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w      a    tangential  velocity  relative  to  the  rotor 

v       =    angle  between  axial  direction  and  absolute  velocities 

p       =    angle  between  axial  direction  and  relative  velocities 

R      3    tip  radius 

u      =    tip  speed 

uu      =    angular  velocity  of  rotor 

r       =    radius  measured  from  compressor  centerline 

p       =    static  pressure 

p       =    density  (assumed  constant  throughout) 
Subscripts  refer  to  stations  as  shown  in  Figure  1. 

In  "two- dimensional"  flow  through  an  axial  compressor,  the 
streamlines  are  assumed  to  lie  on  cylindrical  surfaces  with  constant 
radii.     The  pressure  gradients  in  the  radial  direction  are  ignored,  al- 
though it  is  shown  later  that  at  design  for  vortex  blading,  the  flow  is 
two-dimensional  and  the  pressure  gradients  are  properly  accounted  for. 
Since  an  important  difficulty  in  applying  two-dimensional  flow  theory  to 
this  design  arises,   we  will  discuss  it  briefly  here. 

The  blades  are  of  the  free  vortex  design.     Thus 

c        «    1/r    . 
ul 

A  dimensioiiless  radius  is  defined  by 

£    =    r/R    . 

Therefore,  the  tangential  component  of  velocity  through  the  stafcor  is 

taken  as 

c        =    ua/e     .  (1) 

ul 

where    a    is  an  arbitrary  constant  to  be  determined. 

The  pressure  drop  across  a  stator  row  is  given  by 


Pl  *  Po  1,2  2 
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=     He*-  c   *)    .  (2) 

p  a  l 

Since    r    is  a  constant  along  a  streamline,  the  Bernoulli  equation  for 
rotating  axes  gives,  for  the  pressure  rise  across  a  rotor  row, 
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(w.      -  w.,    )    .  (3) 


p  '1  2 

The  pressurerise  across  a  stage  is  then 

2         o  i*2  2.2  2.  ... 
T-    =    A(w.      -  w,     +c       -  c.    )  .                     (4) 

p  -  v    1  Z  a  1    '  x   ' 

At  design,   c       =  0  ,   and  from  the  velocity  diagrams  we  obtain 
U2 

=    u  «•  +  u£ 
u1  %      "* 

W  =     u£ 

U2 
2  2,        2  2  ,      2  ,a  ,  *  .2 

W ,         =     C         +  W  =C         +  U      ( 7T  +  s ) 

1  a  u.  a  x|      " 

w_      =c       +w        =c       +  u  £ 

2  a  u2  a  ^ 

2  2  2  2         2  .a  .2 

c,       =    c        +  c         =    c        +  u    (tt) 
1  a  u.  a  | 

After  substitution  of  these  relations  into  (4),  we  obtain,  for  the  pres- 
sure rise  across  the  stage, 


(5) 


(6) 

Thus,  when  the  compressor  is  operating  at  the  design  point, 

■■'    =    2a   . 
However,     a    is  not  unrestricted.     It  has  been  found  experimentally 
desirable  (E.ef.    1)  to  have 
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where    p    is  between  0.  5  and  0.  6  .     Oefining  a  non-dimensional  flow 
coefficient  as 

<P    ■    ca/u    .  (8) 

tliis  restriction  can  be  written 

(*;*i2>2    ,„*L  <  1  ,  (9) 

4r/i-{3       iz/p 

This  condition  must  be  satisfied  over  the  entire  length  of  the  blade. 
Since  the  restriction  is  most  severe  at  the  hub,   values  of    $    and    <j> 
evaluated  there  will  be  conservative  over  the  rest  of  the  blade.     An  ex- 
isting experimental  set-up  prescribed  the  hub  ratio    |,  =  0.  6  .     There  is 
another  restriction  on  the  velocities  that  is  not  easily  described  in  analy- 
tical form.     This  is  the  amount  of  turning  that  can  be  accomplished  ef- 
ficiently in  the  blade  rows.     Due  to  lower  blade  speed  at  the  hub  and  the 
adverse  pressure  gradient  in  the  rotor,   this  condition  also  is  most 
severe  at  the  rotor  hub.     p    was  chosen  as  0.  6,  and  values  of  the  turning 
angle  were  obtained  for  different  values  of    '   .     The  quantities    $  =  0.  70 
and    <{>  *  0.  58  were  subsequently  chosen.     The  value  of    #    can  be  raised 
substantially  for  a  compressor  with  hub  ratio  in  the  range  0.  7  to  0.  75. 

So  far,  no  difficulty  with  the  two-dimensional  theory  has  arisen. 
This  is  not  unexpected  since,  as  mentioned  before,  vortex  blading  is 
exactly  two  dimensional  at  the  design  flow  rate.     However,  we  nave  yet 
to  consider  the  off-design  conditions  where  a  radial  component  of  ve- 
locity does  exist.     At  design,  the  work  coefficient  was  constant  over  the 
blade  height.     This  is  desirable  in  order  to  cut  down  losses  due  to  mix- 
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ing  at  the  end  of  the  compressor,   and  actually  is  necessary  on  all  stages 
after  the  first  few  because  a  piling  up  of  pressure  at  some  radius  is 
impossible  for  a  very  large  number  of  stages.     Thus,  this  restriction 
will  be  applied  at  off-design  flow  rates  also,      rhen 

A  1  ""}  **)  O  "^ 

i|i    *  % — m-    ■       2    2   (wl    "w?  +c?   "ci^  3  constant.  (10) 

yp>i-    R  Ui    R 

Remembering  that   c       4  0  for  off- design  conditions,   then 

u2 
w        =    u£  +  c        , 
ul  ul 

wu£    =    ^  +  cu  ,    • 

and fiie  work  coefficient  becomes 

*    =    2  i-  (c      -c      ) 
U        Uj      u2' 

or 

=    24CS  ♦  4>(tai*Y1  -  taap*2)!     .  (11) 

A  first  approximation  can  be  made  that  tanv.  and  tan  6-,  ,  correspond- 
ing to  the  leaving  angles  of  flow,  are  independent  of  the  inlet  angles  and 
hence  are  independent  of  flow  rate.  If  this  is  done,  the  local  value  of  <J> 
at  the  root  is 

»<e  =  o-6>  ■  i.zit°^:LVl)    •  (12) 

Since  under  our  approximation  (  tan0?  -  tanv.)    is  constant  and  is  great- 
er than  zero,  we  see  that  for    *  >  0.  72  the  flow  coefficient  at  the  root 
becomes  negative  and  is  meaningless.     If   tan^-    and    tanv.    are  allowed 
to  change  with  flow  speed  and  this  change  calculated  by  cascade  theory, 
the  singularity  in    6    is  not  relieved  appreciably.     Since  values  of  the 
ideal  pressure  coefficient  greater  than  0.  72  would  be  required  in  starting. 
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for  instance,  the  two- dine  nsional  analysis  of  the  flow  is  unsatisfactory. 
It  must  be  mentioned  that  if    ♦    could  be  pushed  as  high  as  0.  72  in  a 
standard  compressor,  there  would  be  the  same  difficulty.     Normally, 
axial-flow  air  compressors  have  an  ideal  pressure  coefficient  of  be- 
tween 0.  4  and  0.  5  .     It  will  be  shown  that  a  three-dimensional  flow 
analysis  of  the  compressor  does  not  contain  this  singularity  in    $  . 
This  difficulty  is  mathematical  rather  than  physical,  and  appears  to  be 
the  result  of  an  invalid  linearization  of  the  three-dimensional  equations 
of  motion.     There  seems  to  be  no  way  of  calculating  the  ideal  flow  pat- 
terns except  by  a  numerical  procedure  with  a  radial  shift  of  the  stream- 
lines taken  into  account. 

The  assumptions  made  in  the  three-dimensional  flow  analysis 
are: 

(1)  incompressible  perfect  fluid; 

(2)  axisymmetric  flow; 

(3)  equations  are  written  for  conditions  far  upstream  and  far 
downstream  of  a  blade  row. 

The  notation  used  is  the  same  as  in  the  two-dimensional  analysis 
except  that  the  odd  subscripts  refer  to  stations  far  downstream  of  stator 
rows  and  the  even  subscripts  refer  to  stations  far  downstream  of  rotor 
rows  (see  Figure  2). 

The  Bernoulli  equation  for  flow  through  a  stator  row  is 

2p2n+l    ^2                  2                 2Pza  2  2  ..,. 
+  c             +  c               =   +  c         +  c            .                 (13) 

P  U2n*l         a2n+l  P  u2n        a2n 

The  Bernoulli  equation  for  flow  through  a  rotor  row  is 
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The  equations  for  the  pressure  gradient  in  the  radial  direction  are 

2  ,  2 

dp2n    _.        Sn                    .           d?2n+l    _        Sn+l  ,.-. 

^_   =    p  _ and  3- «    p  - .      (15) 

2n  2n  2n+l  2n+l 

By  using  (15)  to  eliminate  the  pressure  terms,   equations  (13)  and  (14) 


become,   respectively, 
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Considering  a  cylindrical  stream  tube,  the  continuity  equation  can  be 

written  in  the  differential  form 

c        r,   dr,      =    c  r-,    .  ,dr~    ,,    . 

a-      2n     2n  a2  +1    2n+*      2n-M 

Within  the  scope  of  the  assumptions  previously  mentioned,  these  equa- 
tions are  quite  general.     It  is  now  assumed  that  the  leaving  angles  of 
both  stator  and  rotor  are  constant  with  respect  to  flow  rate  and  give 
vortex  flow  at  the  design  flow  rate.     Thus,  evaluating  the  leaving  angles 
at  the  design  condition,  indicated  by  subscript    d  ,  we  obtain 
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or,  in  the  non-dimensional  notation, 
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'd  l 

1         *d§2n+l         "^d    *2n+l 

Similarly,  we  obtain  for  the  rotor  leaving  angle 

ta^2  ■  Tr   • 


(18) 


(19) 


(20) 


With  the  leaving  angles  so  specified  and  the  following  definitions, 

a  =  75^       and  b  x   ^d  • 

the  equation  for  the  flow  through  the  stator  becomes 
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The  equation  for  flow  through  a  rotor  row  becomes 
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(22) 


a  +&2n-l 
The  continuity  relations  are 

*2n-lS2n-ld*2n-l    *    *2n^2nd62a   =    ^n+l^n+l^+l  *       (23) 
The  two  pairs  of  non-linear  equations  (21)  and  (22),   and  (23)  can  be 
solved  approximately  for    $       as  a  function  of    £       for    m  =  0,  1,  2, . . . 
in  succession. 

Notice  the  particularly  simple  solution   <j>      =  A,    for  all  values 
of   m  .    It  follows  that  there  is  no  streamline  shift  and  no  radial  velocity 
in  this  case,  and  the  analysis  of  the  flow  reduces  to  the  two-dimensional 
one  previously  considered.     It  is  largely  because  of  this  analytical  sim- 
plification that  the  free  vortex  design  has  been  used  so  frequently.     How- 
ever, for  off-design  flow  rates,  no  such  simple  solution  exists,  even 
for  the  free  vortex  design. 

To  illustrate  the  method  of  solution  to  these  equations,  let  us 
consider  an  off-design  constant  axial  flow   q>      approaching  the  first  sta- 
tor.     In  this  case,   care  must  be  taken  in  applying  equation  (21),   since 
the  assumption  of  a  whirl  velocity  is  incorporated  in  the  right-hand  side 
of  this  relation.     For  this  incoming  flow,  equation  (13)  becomes 

2Pl  ?  ?  2P, 
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Since    dp   /dr    ■  0  ,  equation  (16)  becomes 
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Substituting  from  equation  (19),   the  velocity  diagram  relations,   and  the 
definition  of    a  ,  the  equation  corresponding  to  equation  (21)  is 


2.    dA, 
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(••*>*• 


or 


$.    «    constant. 

This  constant  is  then  determined  by  integrating  the  continuity  relation 

from  the  hub  to  the  tip. 
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Equations  (21)  and  (22)  are  applicable  to  all  of  the  other  blade 


rows  now.     The  relation  for    <j>2    becomes,  using  equation  (22), 

,2.2. 

d 


1 — 


b1 


2^2 


Integrating,  we  obtain 


K  ♦  £  f 

*  2 

^2    =    b~"7 7" 

*     +  « 2 

where    K    is  determined  as  before  by  integrating  the  continuity  expres- 
sion. 

The  flow  coefficients  for  the  following  blade  rows  are  obtained  in 
the  same  manner;  however,   the  integration  for  these  must  be  done  nu- 
merically.    Also,  another  mathematical  difficulty  arises  in  the  solution 
for    $3    and  the  remaining  stages.     The  existence  of  both  independent 
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variables  on  the  right-hand  side  of  both  {21}  and  (22)  will  require  an 
iteration  process.     The  procedure  is  to  first  let    £2      .  *  £-      and  solve 
for    $2     .  .     The  continuity-  relation  can  now  be  integrated  from    £.    to 
g  ,   and  a  relation  between    §2«+l    an<^    6?      obtained.     Substitution  of 
this  relation  into  the  differential  equation  for    4>2    , ,    yields  a  refined 
value  for    4>2n+l  *     Fortunately,  the  correction  to  the  flow  coefficient  is 
not  large,  and  the  iteration  process  is  generally  not  necessary. 
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UI.     APPLICATION  OF  THE  THREE-DIMENSIONAL  THEORY 

The  calculation  of  the  flow  field  as  it  progresses  down  the  com- 
pressor was  carried  out  for  two  off-design  conditions.     The  blades  were 
chosen  such  that    (j>,  ■  0.  57751    and    ♦  ,  m  0.  70000  .     The  average  flow 
rate  for  the  off-design  conditions  was  taken  as    <{>  =  .  46000  and    $  * 
.  69000  .     Figures  3  through  6  show  the  variation  of  the  flow  coefficient 
versus  dimensionless  radius  at  various  stations  in  the  compressor. 
The  flow  field  at  lower  than  design  flow  rate  is  converging  fairly  rapidly 
to  a  steady-state  condition,  while  at  the  higher  flow  rate  the  convergence 
does  not  appear  to  be  as  good. 

The  error  in  choosing   £-  +i  =  Ipn    as  a  **rst  approximation  is 
shown  in  Figure  7.     The  calculation  was  carried  out  for  the  fourth  rotor 
row.     An  error  of  only  approximately  2  percent  of   <j>g    at  the  mid-height 
of  the  blade  is  introduced  by  the  first  approximation.     It  is  felt  that  real 
fluid  effects  will  be  greater  than  this  amount. 

The  local  values  of    *    for  the  two  off-de3ign  flow  rates  are  shown 
in  Figure  8.     Notice  that    '     is  very  nearly  independent  of  the  radius, 
agreeing  with  our  assumptions  leading  to  equation  (10). 

In  standard  compressors  there  exists  a  perturbation  solution  to 
the  three-dimensional  equations.     There  it  is  assumed  that  the  flow   <J>. 
from  a  stator  is  approximately  equal  to  the  flow    <J>_    from  a  rotor.     Thus 
the  procedure  is  to  find  a  solution  for  which    ;^2-<p.  '  «  (<J>2+<l>i)  •     To 
the  zeroth  order,  this  results  in  a  linear  compressor  performance 
curve  in  the  form 

*    s    A  -  B<f>    .  ( 24) 
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The  two-dimensional  theory  also  produced  an  equation  of  this  type, 
since  in  that  case    4>.    was  taken  identically  equal  to    4>.  . 

To  examine  the  possibilities  of  this  approach  for  the  proposed 
helium  compressor,   consider  the  approximate  performance  equation 

*    =    1%Z  +  ^jtanYj  -  |d>2tan|32]    .  (25) 

For  free  vortex  blading,   equation  (25)  becomes 

=    Z    7T-   *1  -t-*?"^2  '  <26) 

26d       1       c|>d     2 

If  we  now  assume 

*x    s    l(*1H2)  +  iC^|-42)    3    *0  +  A*  (27a) 

and 

4>2   =    t<*i+4>2)  -  i<4>1-^z)   a    <>Q  -  A<>    .  (27b) 

then  equation  (26)  can  be  written 

*  •  2e2  +  J[r*d-2623  +  x^d**2i  •  <28> 

Yd  d 

Typical  values  of    *  ,    and    A<|>    for  a  conventional  compressor  allow  the 

term  containing   A<|>    to  be  ignored.     However,  for  the  compressor  pro- 

2 
posed  in  this  report,     if  ,    is  approaching    2£"    at  the  hub,   and  the  term 

containing    Acj>    can  no  longer  be  neglected.     Figure  9  shows  the  per- 
formance curves  of  both  types  of  compressors  with  and  without  the  in- 
clusion of  the    A<j>    term.     The  inadequacy  of  the  two-dimensional 
solution  in  the  case  of  the  helium  compressor  is  obvious.     To  the  scale 
of  the  graph,   no  difference  can  be  observed  for  the  conventional  com- 
pressor. 

Figure  10  shows  the  total  pressure  variation  along  the  blades 
after  the  flow  has  essentially  settled  into  a  steady  pattern.     For  lower 
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than  design  flow  rates,   more  work  is  being  done  by  the  tips,  while  at 
higher  flow  rates,   the  roots  give  the  greater  rise.     The  cumulative  dif- 
ference after  six  stages  is  still  less  than  half  the  design  pressure  rise 
for  a  single  stage. 

It  was  mentioned  previously  that  it  is  desirable  to  have 

2  2 

wl     "  w2 

»i  <    0.  6    . 

wl 
This  quantity  is  plotted  in  Figure  11.     Possible  difficulties  are  indicated 
at  lower  than  design  flow  rates.     A  higher  hub  ratio  would  eliminate  this 
difficulty.     The  corresponding  ratio  for  the  stator  is  plotted  in  Figure 
12.     Since  the  stators  are  in  a  field  of  decreasing  pressure,  no  diffi- 
culties are  expected  here. 

The  off-design  entering  angles  for  the  stator  and  rotor  are  plot- 
ted in  Figures  13  and  14,   respectively.     The  stators  are  very  critical 
and  will  require  care  in  their  design.     Sharp  leading  edges  of  the  sta- 
tors would  probably  cause  high  losses.     The  off-design  flow  incidence 
angle  deviations  for  the  rotors  appear  satisfactory. 
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IV.     CONCLUDING  REMARKS 

It  has  been  demonstrated  that  high  reaction  blading  can  be  used 
in  multi-stage  axial  compressors  to  produce  a  significant  increase  in 
pressure  rise  per  stage  over  conventional  blading.     This  type  of  blading 
could  be  used  in  a  clo3ed-cycle  helium  compressor  and  greatly  reduce 
the  total  number  of  stages  required.     The  analysis  also  indicates  that 
a  hub  ratio  greater  than  0.  6  would  possibly  be  a  more  efficient  design, 
since  stator  incidence  flow  angle  deviations  become  quite  large  near  the 
hub  in  the  example  shown. 

It  was  also  shown  that  the  usual  linearization  of  the  three- 
dimensional  equations  of  flow  is  not  valid  for  this  type  of  blading.    How- 
ever,  a  numerical  solution  allowing  for  the  shift  of  the  streamlines  does 
converge  to  a  steady,   repeating  flow  pattern  after  the  first  five  or  six 
stages  for  flows  as  much  as  twenty  per  cent  off  design. 
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Figure  3.     Stator  Axial  Velocity  Coefficients,    <J>  =  .  46000  . 
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Figure  6.     Rotor  Axial  Velocity  Coefficients,     <J>  =  .  69000  . 
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Figure  8.     Local  Value  of    $     Versus  Non- Dimensional  Radius. 
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Figure  10.     Total  Pressure  Variation  with  Radius. 
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Figure  11.     Kinetic  Energy  Conversion  Through  the  Rotor  Versus 

Radius. 
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Figure  12.      Kinetic  Energy  Conversion  Through  the  Stator  Versus 

Radius . 
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Figure  13.     Stator  Flow  Entering  Angles. 
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Figure  14.     Rotor  Flow  Engering  Angles. 
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Figure  14  (continued).     Rotor  Flow  Entering  Angle6. 
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FART  II.     A  DESIGN  STUDY  OF  A  MULTI-STAGE  AXIAL 
COMPRESSOR  WITH  BLADING  OF  HIGH  ASPECT  RATIO 
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ABSTRACT  FOR  PART  II 

A  method  for  designing  blades  of  high  aspect  ratio  in  a  multi- 
stage axial  compressor  is  investigated.     This  method  requires  a  cal- 
culation of  the  flow  induced  by  all  blade  rows.     The  method  is  illustrated 
for  the  limiting  condition  of  hub  ratio  equal  to  unity  and  then  is  carried 
out  for  the  more  realistic  case  of  finite  hub  ratio. 

An  example  of  a  blade  design  for  a  particular  flow  is  carried 
out.     The  results  are  compared  with  those  of  a  previous  design  for  the 
same  flow  conditions  but  based  on  a  theory  applicable  for  blading  of  low 
aspect  ratio. 
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SUMMARY  FOR  PART  II 

The  flow  through  moat  present-day  axial  flow  compressors  is 
frequently  analyzed  by  applying  a  three-dimensional  theory  that  is  based 
on  the  assumption  of  blading  with  low  aspect  ratio.     Even  though  existing 
compressors  do  not  have  what  could  be  called  low  aspect  ratio  blades, 
the  theory  predictions  and  experimental  results  show  reasonable  agree- 
ment. 

Where  weight  considerations  are  of  primary  concern,  the  possi- 
bility of  reducing  the  size  of  a  compressor  by  shortening  its  length  is 
attractive.     There  is  no  known  reason  why  the  blade  rows  cannot  be  held 
to  their  same  height  but  shortened  in  chord.     There  is,   however,   a  dif- 
ficulty in  applying  the  aforementioned  analysis.     There,  it  is  assumed 
that  equilibrium  flow  patterns  exist  at  the  leading  and  trailing  edges  of 
the  blades.     THs  is  identical  to  assuming  low  aspect  ratio  for  the  blades, 
since  they  must  be  wide  enough  to  allow  all  changes  in  the  flow  to  occur 
entirely  within  the  blade  row  itself.     Several  extensions  of  the  theory 
have  been  developed  to  take  account  of  the  non-equilibrium  flow  near  the 
blade  rows.     Some  of  these  analyses  assumed  small  vorticity  throughout 
the  flow  field  and  also  a  small  change  in  vorticity  through  the  blade  rows. 
Another  theory  relaxed  the  restriction  of  small  vorticity  throughout  the 
flow  field. 

The  theory  presented  in  this  thesis  is  an  extension  of  the  latter 
approach,   and  is  based  on  the  assumptions  of  closely  spaced  blades  of 
high  aspect  ratio.     The  equations  of  motion  are  linearized  by  a  perturba- 
tion scheme  that  assumes  small  changes  in  perturbation  velocities  across 
the  blade  rows.     To  illustrate  the  method,  the  calculation  is  first 
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restricted  to  the  limiting  case  of  hub  ratio  close  to  unity.     Then  the 
more  involved  problem  for  any  hub  ratio. is  investigated.     An  example 
is  worked  for  blades  of  aspect  ratio  about  three.     Here,  there  appeared 
little  difference  between  the  theory  applicable  to  blades  of  low  aspect 
ratio  and  the  theory  of  blades  of  low  aspect  ratio. 
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I.     INTRODUCTION 

The  theory  of  three-dimensional  flow  in  the  axial  compressor 
described  in  Part  I  of  this  thesis  is  concerned  only  with  flow  conditions 
far  upstream  and  far  downstream  of  a  blade  row.     This  theory  would  be 
applcable  to  blades  of  very  low  aspect  ratio.     The  difference  between 
these  velocity  profiles  and  the  actual  profiles  near  the  leading  and  trail- 
ing edges  of  the  blades  was  assumed  to  be  negligible.     It  was  on  this 
basis  that  the  entering  flow  angles  for  the  blades  were  determined.     If 
the  aspect  ratio  is  not  very  low,   it  is  possible  to  obtain  more  exact  in- 
formation about  the  flow  pattern  near  the  blade  rows  by  making  use  of 
some  of  Marble's  results  (Ref.    1).     This  has  been  done  by  Bowen, 
Saber  sky,   and  Rannie  in  Ref.   2.     The  method  carried  out  by  these  in- 
vestigators forms  the  basis  of  the  theory  to  be  discussed  in  this  report. 
The  approach  is  to  consider  the  disturbances  due  to  rows  of  blades  of 
high  aspect  ratio  placed  very  close  together.     Here,   the  flow  in  any  par- 
ticular blade  row  is  influenced  by  neighboring  blade  rows.     In  the  three- 
dimensional  theory  of  Part  I,  the  changes  in  the  flow  pattern  were  as- 
sumed to  occur  entirely  within  the  blade  row  itself  and  one  blade  row  did 
not  affect  its  neighbor. 

The  analysis  begins  with  the  calculation  of  the  disturbances  due 
to  an  infinite simally  thin  disk.     These  disturbances  can  then  be  super- 
imposed to  form  a  blade  with  finite  axial  extent.     The  effect  of  many 
blade  rows  can  be  determined  by  the  superposition  of  the  disturbances 
due  to  the  finite  row.     In  this  way,  the  flow  pattern  for  a  multi-stage 

axial  turbomachine  can  be  constructed.     For  a  repeating  flow  pattern, 
the  summation  of  the  effects  of  all  blade  rows  is  possible. 
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n.     TWO- DIMENSIONAL  THEORY 

A.     Equations  and  Boundary  Conditions 

Before  proceeding  directly  into  the  three-dimensional  theory  of 
high  aspect  ratio  blading,  it  ie  informative  to  investigate  first  the  limit* 
ing  case  of  high  hub  ratio.     This  problem  is  simpler  because  the  radial 
pressure  gradients  vanish  far  from  the  blade.     The  analysis  is  based 
upon  the  disturbances  due  to  a  disk  at    x  a  0    (see  Figure  l)f  and  the 
equations  of  motion  are  written  subject  to  several  assumed  conditions. 
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Figure  1 


These  assumptions  are:    (1)  there  are  no  circumferential  variations; 
that  is,    •*—  (   )  a  0  •     This  is  equivalent  to  having  a  disk  with  an  infinite 
number  of  blades.     (2)    The  flow  is  incompressible,  the  density    p    being 
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held  constant.     The  continuity  and  momentum  equations  then  become 

du    ,    3v 


"§x  +  "Sy 


«    0  (1) 


Du    .       8u  1   Dp  .,. 

U3£+V3y    =    "35  <2) 

3v    .        3v  1   dp  /al 

U  "53E  +  v  "5y    s   °  <4) 

where    u  ,  v  ,  and   w    are  the  gas  velocity  components  in  the    x  ,  y  ,  and 
z    directions,  and   p    is  the  static  pressure. 

The  boundary  conditions  to  be  imposed  on  this  system  are  as 
follows: 

(1)  v=0    at   y=0    and   y  a  b    for  all    x  ,   since  the  walls  of  the 
compressor  are  streamlines. 

(2)  u   is  continuous  at    x  =  0    from  continuity  through  the  actu- 
ating disk. 

(3)  v    is  continuous  at    x  =  0    due  to  the  component  of  force 
being  zero  in  the    y    direction. 

(4)  w    is  allowed  to  change  dis continuously  across    x  =  0  . 

(5)  At    x  =  -  co  , 

u   =   u(y)  -  cuo(y), 

v   a    0  , 

w    m    w(y)  -  rwQ(y)  . 

p    =    p  -  e.p      *    constant. 

(6)  At  x  a  +  oo  , 

u    a    u(y)  +  ^uo(y)  , 
v    a    0  , 
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w    *    w(y)  +  ewQ(y)  . 
p    *    p  +  cp      =    constant. 

These  last  two  boundary  conditions  are  written  in  a  form  so  as  to  con- 
form with  the  following  perturbation  scheme.     Far  upstream  and  far 
downstream  of  the  disk,  the  flow  is  independent  of   x  ;  any  disturbance 
due  to  the  disk  having  been  smoothed  out.     Close  to  the  disk,  the  flow  is 
disturbed  by  a  small  amount,  and  the  following  perturbations  are  as- 
sumed, where  the  symbol    e    is  carried  only  to  denote  the  assumed 
small  quantities. 

u   m   u(y)+    cuQ(y)  +  cujfx.y) 

v   =    e  Vjfoy) 

w    a    w(y)  +  cwo(y)  +  cw1(x,y) 
p    *    P  +  ?P0  +  epjCx.y) 

Notice  that    p    and   p      are  assumed  independent  of  either    x    or    y    and 

that  the  mean  value  of   v    is  zero.     The  top  signs  are  for    -  co  <  x  <  0 

and  the  bottom  signs  are  for    0  <  x  <  +  co  . 

Introducing  equation  (5)  into  the  equation  of  motion  yields  the 

linear  set 

8u.         9v. 

-35T  +  ST    '   °  <6> 

_  8ui       8u  „  1  *»l  m 

u  "53T  +  3y  vl   =   "  p  1*  <7) 

UTZ-    *   *7"Sy  (  ' 

9w  — 

«^r+-Hvi  =  °  <9> 
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B.     Solution  of  the  Equations  for  the  Actuating  Disk 

The  first  three  of  the  preceding  equations  (6,   7,   and  8)  can  be 
combined  to  give  the  differential  equation  for    v  ,   as 


a2v,       92v,       ,  n2_ 

1     ,  1        1    0  u  n 

TT  +  "T~T  "~7^  vl    =    ° 
dx  dy  u  8y 


(10) 


with  the  boundary  conditions  that    v  =  0    at   y  =  0    and   y  =  b  .     A  solu- 
tion to  equation  (10)  that  satisfies  these  boundary  conditions  is 


co 


TT     V     A  1  ^    b" 

'l    =    U   L  An  e 


V*> 


n=l 


where 


(ID 


U    » 


u  dy 


and    Y  (y)    satisfies  the  differential  equation 


d2V 


n 


dy' 


/  Xn  \       U7' 

Vr")  "  -=■ 


Y      3    0 
n 


with  boundary  conditions 


Y  (o)    ■    Y  (b)    =    0    . 
n*   '  n*   ' 


Equation  (6)  then  yields 


«,      =     + 


co 

UEAn 
n=l 


x 


ie±X^Yn(y) 
Xn  n 


(12) 


where  the  prime  denotes  differentiation  with  respect  to   y  .     Notice  that 
the  boundary  conditions  on   u    are  satisfied  at    x  =  -  oo    and   x  =  +  co  . 
?rom  the  continuity  of   u    across  the  disk,  the  following  relation  is 
obtained 


oo 


uo(y)   .    -  U  £  Aa  ±  yyy)  . 


(13) 


n=l 
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Since    u    and   u      are  to  be  chosen  and  the  shape  of  the  blades  is  to  be 
determined,   it  is  possible  at  this  point  to  choose    A    s  0    for    n  >  1  . 
Also,  if   u   is  chosen  as  a  linear  function  of   y  ,  then 

Y,    =    «n-5- 


and 


Xj    c    w 


Note  that  it  is  not  necessary  that  these  particular  choices  be  made.     It 
is  done  here  only  to  reduce  the  amount  of  algebra  in  the  rest  of  the  analy- 
sis.    It  will  be  pointed  out  later  that  these  choices  can  indeed  represent 
the  flow  through  a  conventional  compressor  quite  well.     Therefore,  the 

perturbation  velocities  become 

x 

J-     * 

u,    s    +  u     e 

1        —     o 
and 

v,    ■    -  e  I  u  (n)  dn     . 

*  Jo 

0 

Using  equation  (9)  »     w.    is  determined  as 

Wl     "    l-r^e"*^      I     U(T))dTl     ,  (14) 

u  0 


and  with  the  use  of  equation  (7), 


1  ,    ±wb- 

p*l    =    +e 


y 

—  du  i    \  i 

u  u    -  -5—       u  (n)  dv 
o      ay    J     ox 


(15) 


0 

There  are  several  requirements  that  we  can  now  impose.     First, 
we  can  require  the  total  pressure  rise  across  the  disk  to  be  a  constant. 
This  will  insure  that  the  flow  pattern  will  repeat  itself.     This  is  a  similar 
condition  to  that  set  upon  the  work  coefficient  in  Part  I  of  this  paper.     To 
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calculate  the  change  in  total  pressure,   one  needs  the  change  In  whirl 
and  the  change  in  static  pressure.     The  change  in  wnirl  is 

Aw    a    w(+o,  y)  -  w(-o,  y)  , 
and  with  the  use  of  equations  (5)  and  (14)  becomes 

Aw    a     Zcw     -  2c  —  -3 —     i   u  (r)  dv.    . 

u        o 

The  change  in  static  pressure  is 

^    *   ?  <+o.y)-£  <-o,y)    . 
and  using  equations  (5)  and  (15),  becomes 


£&    =    2,  £2  +  2c 
P  P 


u  u    -  -t^      u  (n)dn 

o       ay    J      <r  "     ' 
0 


Since    u    and   v   are  continuous  through  the  disk,  the  change  in  total 

pressure  is 

ApT         Ap       - 
L    =  .££  +  w  Aw 

P  P 

or 

_  »  2c  -.  +  2c  (u  u0  +  w  wo)  +  2c  (^  +  -  ^  )  j  uo(n)dTi    ,        (16) 

u  0 

and  we  require  this  to  be  constant.     Thus,  we  now  have  the  four  unknown 

functions  of   y  :    u  ,  u    ,  w  ,  and    w    ,  and  two  constants    p    and   p 

o  o  o 

with  one  equation  (16)  relating  them.     Also,   since  we  are  interested  only 
in  the  pressure  rise  in  a  compressor,  v/e  are  not  immensely  concerned 
with  the  mean  pressure  ,  p  .     The  pressure  change  across  the  disk  is 
2p    ,  and  this  is  the  more  important  quantity. 

Another  condition  that  v/e  can  require  is  that  the  total  pressure 
far  upstream  and  far  downstream  be  constant  from  root  to  tip.     Note 
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that  this  is  not  a  necessary  condition,  but  is  a  common  choice  in  com- 
pressor design  and  gives  two  more  relations  connecting  our  four  un- 
known functions  of   y  .     For  the  total  pressure  to  be  constant  we  must 
require  that 


,2      , - 
I     +  (w- 

and  (17) 


TJ  —  2  —  Z 

*-  +  i  [(u-pu  )    +  (w-ew  )         =    constant 
p        A    '  o  o 


2  2 

*•  +   aC(u+"u  )    +  (w+pw  )    ]    ■    constant 
p  o  o 

and  since    p/p    is  constant  far  upstream  and  far  downstream  of  the 

disk,  we  derive  the  following  two  equations  by  adding  and  subtracting 


equations  (17)  and  keeping  only  first  order  in    c  : 

(18) 


U       +  W         a     K. 


U    U  +  W    W     a     K_ 

o  o  c 

where    K.    and    K,    are  undetermined  constants. 

The  second  condition  of  equations  (18)  is  now  actually  the  same 
as  equation  (16).     If  the  disk  were  to  act  as  a  stator,  then  (16)  would 
reduce  to 

/du    ,    w    dw  \     |\ 


-  T  *  uuo  *  wow  -  Vdy  +  r  -ay  /  J  %{v)dr]  '  (19> 

u  0 


Since    p/p    is  a  constant,   and  from  (18)    uu    +  w   w  a  K,  ,  the  last 
o  r  o         o  Z 

term  of  equation  (19)  must  be  a  constant.     Since  the  integral  is  obviously 
zero  at   ysO,  then  this  constant  must  be  zero.     Thus,  with  the  assump- 
tion of  constant  total  pressure,  we  have  four  unknown  functions  of   y 
and  two  conditions  upon  them.     Hence,  for  example,  we  could  pick   u  , 
u    ,     and   p    ,   and  the  flow  will  be  completely  determined.     Once  this 
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has  been  done,  the  chape  of  the  blades  can  be  found. 

To  derive  the  effect  on  the  flow  due  to  a  blade  of  finite  width    a 
and  height   b  ,  we  can  sum  from    x  =  0    to   x  a  a   the  disturbance  due  to 
a  disk  located  at   x  =  §    whore    £    ranges  from    0    to    a  .     Notice,  how- 
ever, that  here  one  must  make  a  choice  of  the  shape  of  the  loading  on 
the  blade.     That  is,  we  must  require  that  each  disk  contribute  some 
fraction   f(£)d£    of  the  disturbance  of  the  blade.     In  the  problem  pro- 
posed in  this  report,  we  are  at  liberty  to  pick  the  function   f(|)  .     We 
then  go  on  to  And  the  shape  of  the  blade.     However,  the  inverse  problem 
of  specifying  the  blading  and  then  finding  the  flow  is  much  more  diffi- 
cult.    In  that  case,  the  loading  function  is  in  general  a  function  of  both 
x   and   y    and  depends  on  the  entering  flow  angles. 

C.     Superposition  of  Actuating  Disks  for  Multi-Stage  Compressors 

Given  the  function   f{|)  ,  we  can  proceed  to  determine  the  ve- 
locities   u    and   w   throughout  the  compressor.     The  function   f(|)    is 

normalized  such  that 

a 

r«t«  ■  i. 

V 

0 

For  a  finite  stator  of  width    a  ,  height   b  ,  and  leading  edge  at   x  «  0  , 

the  axial  velocity  for    x  >  0    becomes 

.      x    a  x 

usu-eu+eu     e  f(|)  e  d§    .  (20) 


for    0  <  x  <  a  , 

-            \  P 

Co 

"i./*w 

a 
<1£  -    ff(fc) 

o 
X 

ir 


1  -  e     ° 


d£(  •     (2D 
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For   x  >  a  , 

x     a  * 

u   a    u  +  cuQ  -  GUQ  e       ^  J  f(£)  e       ^  d£    .  <22) 

0 
The  equations  for  a  rotor  are  similar  with  only  the  algebraic  sign  of 
u      changed  throughout. 

Equations  (20)  -  (22)  give  the  flow  due  to  one  blade  row  alone. 
For  a  repeating  flow  pattern,  the  velocities  are  the  same  as  if  there 
were  an  infinite  number  of  blade  rows.     Thus,  for  the  complete  flow  pat- 
tern, we  must  superimpose  the  flows  due  to  many  stators  and  many  ro- 
tors.    The  following  definitions  will  be  used: 

a„    *   width  of  a  stator  row 

8 

a      a    width  of  a  rotor  row 

r 

c      «    gap  between  stator  and  rotor 

c      =    gap  between  rotor  and  stator 

f  (£)dg    a   fraction  of  disturbance  due  to  stator 
s 

f  (£)d£    s    fraction  of  disturbance  due  to  rotor 

<t    =    a    +  c    +a    +c 
s        s        r        r 

The  leading  edge  of  the  stators  are  then  set  at 

x   *    ml         where        m  *  0,   +1,  +2,   ... 

and  the  leading  edgesof  the  rotors  are  at 

x    ■    a     +  c     +  ml    . 

s        s 

The  axial  velocity  in  a  stator  row  due  to  all  blades  becomes 


X 

p 


u  •  u(y)  +  cuQ(y)  j   j  £Q(i) 

0 
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, .  /tw' 


di  -    ye) 


X 


1  -  e 


+£<*-£) 


d£ 


CD 


-  IT 


x-a  -c  -m£] 


-  e 


^  u0(y)  J*r(t>«    "Fd£  J 


m»0 


co 


a  » 

or  +  irg  -£  fx-nvt] 


■f  €     )     u  (y)        f  (£)  e      M  d|  e 


o 
m=l  0 

CD 


+  c    £    u0(y)   Jfs(|)  e       ^d£  e 
m=l  0 


+£  [x-nvt] 


CO 


-  e 


«o<y> 


a.  «. 


mal 


fs<e)  e       u  d£  e 


Carrying  out  the  summations,  equation  (23)  becomes 

a 


x 


4y)  +  euo(y)j  j   f(|) 
0 


1  -  e 


-?<*-£) 


r 


j 


(23) 


d£-   j  f(£)^l  -e 


+  £(*-£) 


=£ 


x 


™o(y) 


£  (I)  e       *  d£  e    F  e   F     8     s 


1  -  e 


r  +*&        -**    -£(a  +c   )       Ps  -irf 

£r(l)  e       ^dle       beDrr-|  yi)  e       Dd£  e 


x       W  , 


0 

x         w 


a 


(24) 


For  the  axial  velocity  in  a  rotor  row,  change  the  sign  of   u      and  inter- 
change the    s    and    r    subscripts  in  all  terms  of  equation  (24). 
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The  corresponding  calculation  of  the  whirl  velocity  in  a  stater 
results  in 


w  =  w 


<y)  +  e  Jys: 


0 

a 


X 


i  dw  p 


wo<y)  -  r  ^  j  V1^  e 


■5^ 


d£ 


wc(y)  - 


-  W  J  uo(r)dT1  e 


d| 


-    7 

1    dw     P       .    „  1 

-  c  -r  3v     l  n0<T1,dT1 

u  «  ' 

0  l-e 


irT    /  J  fr(&)e  **«        e 

V   (  o 


a  A  a 

£r(ft)e      %e      *e   *     r     r 


-  j  fs(|)  e         d£  e       e 

0 


a                        a  » 

s                   £  x  I 

+    I  £  (|)  e       bd|  e       b  e       b 

J     s 


(25) 


The  whirl  velocity  in  a  rotor  is  determined  by  changing  the  sign  of   u 
and   w      and  interchanging  the    s    and    r    subscripts  in  equation  (25). 

An  example  will  be  shown  later  for  which  the  following  assump- 
tions are  made: 

a     a    a     a    a   . 
s  r  ' 

c      a    c      a    c    , 
s  r 

ys)  =  ye)  =  i/a  , 

t    =    2(a  +  c)    . 


D.     Flow  Patterns  for  Blading  of  High  Aspect  Ratio 

If  then  the  exponentials  arc  expanded  in  polynomials  of    a/b  , 
c/b  ,  or    x/b  ,  as  the  case  may  be,  and  terms  through  second  order  are 


(26) 


x       .       rrx  ,  tta. 


u    =    u  -  r\x  it  (  £  -  -«-   ) 


(27) 
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retained,   the  velocities  reduce  to 

u  =  u  +  c  uo*  (  £  -  £  ; 

w    =    w  +  ew     (2--  l).rr(2i.  1.^  +  ^  ) 

o  *    a        '  *     a  o        b 

for  a  stator,   and 

w    =:    w-cw     (2~-l)  +  Cr(2--l-22E+Z£) 
o  *     a         '  x     a  b        Zb  ' 

for  a  rotor.     In  both  cases 

r  K  -r  t£   I  V^  dl1    •  (28) 

o 

The  entering  and  leaving  flow  angles  for  both  stator  and  rotor 
blades  can  now  be  calculated.     The  defining  equations  for  these  angles 
are, in  the  case  of  a  stator, 

u(x  =  0) 

tan  v ,     =    — t rrr 

'  1         w(x  =  0) 

and  (29) 

u  (x  =  a) 

ten.  Vt    =    — )  ( 

*  2         w  (x  =  a) 

where    y,    is  the  stator  entering  angles  and    v2    is  the  stator  leaving 

angle.     For  the  rotor, 

tan  6      =  u(*-0) 

pl  blade  speed  -  w(x=07" 


and 


t       a  u(x  g  a) 

'  2         blade  speed  -  w(x=a) 


In  summary,   given  the  blade  loading,   blade  aspect  ratio,   the 
axial  velocity  through  the  compressor,   the  whirl  velocity  at  one  radius 
of  a  blade  row,   the  desired  pressure  rise  across  a  stage,   and  the  de- 
sired degree  of  reaction,    one  can  compute  the  shape  of  the  blades  subject 
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to  the  restriction  that  they  be  of  high  aspect  ratio  and  very  closely 
spaced. 
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III.     THEORY  FOR  FINITE  HUB  RADIUS 

The  analysis  for  the  limiting  case  of  hub  ratio  equal  to  one  is 
extended  to  the  axisymmetric  case  of  finite  hub  radius  with  only  slight 
modifications.     The    y    coordinate  goes  over  into  the  radius    r  .     The 
velocity  in  the    r    direction  remains    v  .     The  major  change  here  is  that 
the  radial  pressure  gradients  must  be  accounted  for  properly. 

The  equations  of  motion  in  the  axisymmetric  case  are 

37+7  S?  (rv)    =    ° 

8u    .        8u  1    dp 

uSx-+vU7    =    "IS 

2  <31) 

8v  ,        dv       w  1    8p 

u  ^E  +  v  "57  *  —    s    '  p  « 

9w  8w        vw  n 

u  -k—  +  v  -»—  +  =    0 

Ux  Or         r 

As  before,   we  assume  a  perturbation  approach  of  the  form 

u    =    u(r)  T  cu  (r)  +  eu.(x,  r) 

v    ■    -v.(x,  r) 

(32) 
w  s    w(r)  +  <?w   (r)  +  rWj(x,  r) 

p    =    p(r)  +  cpQ(r)  +  epjfx,  r) 
Notice  that   p    and    p      are  now  functions  of  the  radial  variable,    r  . 


re 

suiting  linear 

equations  are 

9u. 
"5x 

r 

£  (rVl>    =    ° 

u 

8u. 
T5T 

+  v] 

3u 
I  "3T 

i  8?i 

p  "5x~ 

u 

avi 

"3x~ 

-  2 

WW. 

r 

I  8*>l 
P  "BT" 

.„_ 

8w. 

+  v 

8w 
1    "7JT 

v,w 

+  -=—     =    0 

(33) 
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and  also  the  radial  pressure  equilibrium  equations 

2 


and 


}H=^  <«> 


op  2w  v 


1"f  *~W     W 

o  o 

—         lg|i:       '>■■■  S            m«...i«iii>i- 

p   or  r 


which  must  be  satisfied  at    x  -  +  co  . 

Eliminating  the  other  dependent  variables,   the  equation  for    v. 

becomes 

2  2 

a  v.      9  v.      j  dv 

— T  +  -F"1  +  7  ST  "  [-T  +  F<r»i    ■   °  -  <35> 

3x  r 

where 

d    / 1   du\        2     w      d   ,   —  4 

dT  Vr  Hr  /  "  T  "=1 
u  r       u 


_.    .  r     d    /l   du\        2     w      d    .   — v 

F(r)   =  3  dT  (r  3?)  "  "T  =Z  d?  <rw> 


An  appropriate  solution  is  of  the  form 

Vj    =    e        D  R(r)    . 

Substituting  into  equation  (35)  yields 

1  \    2        1 

R"+±R«  +  r(£)     --^-F(r)]R    *    0  (36) 

r 
where  the  prime  indicates  differentiation  with  respect  to    r  .     If    F(r)    is 
chosen  identically  equal  to  zero,   then  (36)  reduces  to  Bessel's  equation. 
If  we  then  take  only  the  first  term  of  the  solution  for    v.  ,   the  boundary 
condition  that    v  =  0    at  the  hub  and  the  tip  gives,   for  the  first  eigen- 
value 

Xj    ■    3.  18 
for  a  hub  ratio  of  0.  6  .     If    F(r)    is  small  compared  to    (X,/b)  '  ,  then 
the  first  eigenvalue  does  not  change  substantially.     We  will  now  choose 
F(r)    identically  equal  to  zero  and  later  compute  its  actual  value  to 
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compare  with    (\./b)     . 

Proceeding  as  in  the  case  with  infinite  hub  radius,  the  perturba- 
tion quantities  become 


x 


u,    =    +  u     e 


±*o- 


+  ^  -     r 

ix    ±xb    r 


I  f\  u0(^)  d*l 


+  x£    r 


w 


1    -    ±3  t-ar  +"7)7  e  .!   ^Vr) 


d-p 


(37) 


p2 
p 


ciu 


r 


=    +  e  [uu0-73?    j   nuo(r)dr) 


where    r.    is  the  hub  radius. 

The  condition  that  the  total  pressure  rise  across  a  disk  be  a 
constant  yields 


Ap, 


P. 


=    2c  —  +  2«:(uu    +ww   )  -  2c 
p  p  o  o' 


1  du" 


dw 


i.    u  j  w   /   w  .  w  \ 
r  clr"      —    VTTF      r"' 


ur 


T}uo(Ti)dn.(38) 


i 


Again  choosing  the  special  condition  of  constant  total  pressure  in  the 
radial  direction  implies  to  first  order  in    e    that 
p  +  1  p  (w2  +  u2)    =    Kx 

(39) 


and 


P„  +  P  (ww     +  uu   )    »    K?  . 


o      r  '       o         _o'         _>2 
Then,   as  in  the  simpler  case,  equation  (38)  applied  across  a  stator  re- 
quires   K2  =  0  .     Combining  equations  (34)  and  (39)  results  in  the  two 
conditions 
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Id.—.  du 

-r  ^F  (rw)    =  -  -ar 

r 

and  (40) 

d    .    2  —       .  2d,—      . 

^  (r   wwQ)    «   -  r    ^  <uu0)    . 

The  remainder  of  the  analysis,   consisting  of  the  summation  of 

the  disturbances  due  tc  the  disks  in  order  to  produce  a  finite  blade  and 

the  summation  over  all  blade  rows  to  derive  the  complete  flow  pattern 

for  the  compressor,  proceeds  exactly  as  in  the  limiting  case,   and  only 

results  will  be  given  here.     With  the  definition 

r 

r   ■   —  (lTF  +  -r)  J  ^o^^    ' 
ru  d 

ri 
the  velocity  components  for  the  stators  are 

U     =     U  +  €Uo\    (^--2b    ) 

and  (41) 

w    =    w+Cwo(2|-l)-,r(2|-l->^+^), 

and  for  the  rotors 

u    =    u-  euo\(£-^) 
and  (42) 

w    a    w-cw0(2£.l)  +  d*(2£-l-^*£). 

The  blade  angles  are  determined  exactly  as  before  with  the  use  of  equa- 
tions (28)  and  (29). 
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IV.     AN  APPLICATION  OF  THE  THEORY 

Before  proceeding  with  an  example  that  will  allow  us  to  verify 
this  theory  against  an  existing  compressor,  we  will  indicate  the  steps 
in  order. 

1.  Choose  the  blade  loading  function  f(£)    for  both  the  rotor  and 
the  stator. 

2.  Choose  the  blade  dimensions. 

3.  Solve  equation  (36)  for  the  eigenf unctions  and  determine  the 
eigenvalues  by  application  of  the  boundary  conditions  on  the  radial  ve- 
locity  v  .     For  many  applications,     F(r)    can  be  taken  equal  to  zero  and 
all  constants  in  the  series  solution  for    v    taken  equal  to  zero  except  the 
first  ones.     For  such  applications,   the  magnitude  of    \.    will  be  approxi- 
mately   ir  . 

4.  Continue  the  summation  procedures  that  result  in  relations 
similar  to  equations  (41)  and  (42). 

5.  Pick   u   and   u    .     Notice    u      must  satisfy  the  condition  that 

o     r  o  7 

o 

,   r  u    dr    a    0 
J         ° 

ri 

where    r.    and    r^    are  the  inner  and  outer  boundaries  of  the  cylindrical 

annulus. 

6.  Determine  w  from  the  first  of  equations  (40).  There  will 
remain  an  undetermined  constant  to  choose. 

7.  Choose  the  pressure  rise  across  a  stage  and  the  degree  of  re« 
action  desired.     The  pressure  rise  across  a  stator    2p      can  then  be  de- 
termined. 
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8.  Determine    w      from  the  second  of  equations  (40).     The  unde- 
termined constant  is  then  evaluated  by  use  of  the  second  of  equations  (39). 

9.  Evaluate  the  function   4>  . 

10.  Calculate  the  axial  and  whirl  velocities  for  both  the  stator 
and  rotor  by  using  the  equations  derived  in  step  4. 

11.  Compute  the  blade  angles  by  use  of  equations  (23)  and  (29). 

It  was  decided  to  use  this  theory  to  compare  with  the  blade  design 
given  in  Ref.  1  for  solid  body  blading.  The  steps  will  be  taken  as  listed 
above  and  the  results  presented  here  and  in  Table  1  on  page  58. 

1.  The  loading  function  was  taken  to  be  constant  for  ease  in 
computation.     A  triangular  distribution  would  be  more  nearly  correct. 

2.  The  ratio    a/b    was  taken  to  be    1/3  . 

3.  F(r)    was  taken  equal  to  zero  and    \.    found  to  be  approximately 

3.  18  .     The  function    F(r)    was  subsequently  checked  and  found  negli- 

.     2 
gible  compared  to    (\/b)    . 

4.  Equations  (41)  and  (42)  were  then  correct  as  stated  in  the 
theory. 

5.  u   and   u      were  then  chosen  to  approximate  as  closely  as  pos- 
sible the  axial  flow  used  in  the  design  of  the  solid  body  blading  in  Ref.    1. 
The  relations  were 

u/uT    ■    -  .  560  |  +  .916 

and 

u  /u_    =    .  185  |  -  .  151 
o      l 

where    £    is  a  non-dimensional  radius  based  on  the  tip  radius  and   u_    is 
the  rotor  tip  speed. 

6.  Equation  (40)  was  used  to  determine    w    in  the  form 
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2 


c 


(i)      =    .342|-.15762--»    . 

The  constant    c,    was  determined  by  choosing    w  ■  .  385    at    £  ■  0.  8  . 

7.  The  ideal  pressure  coefficient  for  the  compressor  was  0.  40. 
Then  by  definition 

ps  .  Pr         0.  40  „    2 
P         P  ^^     T 

where    nm    is  the  blade  tip  speed  that  has  been  used  throughout  this 
analysis  in  order  to  non-dimensionalize  the  velocity  components.     Fifty 
per  cent  reaction  in  this  compressor  occurred  at   £  =  0.  8  .     Since 
2p    =  p    »  we  obtained  the  following  relation 
pQ/p    =    0.  05  uT2    . 

8.  The  equation  for    w      then  became 

Is.   =_L_/052|2-.085|.1°02). 
UT         w/u,/  IF' 

9.  The  function   §   was  evaluated  numerically  and  is  listed  in 
Table  1. 

10.  and  11.     The  values  of   u   and   w   are  listed  in  Table  1  for 
different  values  of   £  .     Similarly,  the  entering  and  leaving  angles    y,  , 
v2  ,  p.  ,   and   p-    are  tabulated.     The  values  of  these  angles  used  in  the 
original  design  are  also  tabulated. 
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V.     CONCLUDING  REMARKS 

For  most  compressor  designs  occurring  in  practice,  the  aspect 
ratio  of  the  blades  is  neither  very  low  nor  very  high.     It  falls  some- 
where in  between.     The  theory  of  Part  I  of  this  thesis  is  applicable  for 
blading  of  low  aspect  ratio,  where  it  is  assumed  that  equilibrium  flow 
patterns  exist  at  the  leading  and  trailing  edges  of  the  blade  rows.     A 
blade  row  is  not  affected  by  neighboring  blade  rows  except  through  the 
incident  velocity  distributions.     The  theory  of  Part  II  for  blades  of  high 
aspect  ratio  contains  the  influence  of  disturbances  of  potential  type  that 
can  be  felt  upstream  as  well  as  downstream. 

The  close  agreement  between  the  results  of  the  two  theories  for 
an  aspect  ratio  of  three  indicates  that  there  the  flow  between  blades  is 
very  near  its  equilibrium  state  at  the  leading  and  trailing  edges. 
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